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Elastic Problem for a Ring of Uniform Force in an 

Infinite Body 1 

William H. Pell 

Kelvin has given an integral representation for the displacements in an infinite elastic 
body produced by a body force acting in an arbitrary portion of this body. By a limit pro- 
cedure, solutions have been found for point and line singularities acting in the interior of 
such a body. A similar method is used in this paper to obtain the displacements produced 
by a uniform force applied along, and normal to, a circle lying in the interior of an infinite 
elastic body. 

1. Introduction 

One of tin 4 important results of the classicial theory of elasticity is the solution of Kelvin 
[1, 3] for the displacements in an infinite elastic solid due to a body force distribution throughout 
all, or a part of the body. From this, one obtains by a limiting procedure the displacements 
for the case in which the body force reduces to a concentrated force at some point in the infinite 4 
solid. This solution is singular at the point of application of the concentrated force, and plays 
a role for the Navier equations of elasticity analogous to that which the potential of a point 
charge plays for Laplace's equation. By the process of superposition, the displacements, due 
to any finite number of concentrated forces with different directions and at different points, 
can be obtained. The principle of superposition and a limiting process such as is used in me- 
chanics to obtain a concentrated moment can be employed to obtain solutions for higher-order 
singularities. For example, one may superpose the solutions for point forces P, a distance h 
apart, and with opposite senses. By allowing h-^0 while Ph remains constant, the solution 
corresponding to a "double force" is obtained. Superposition of three of these with mutually 
perpendicular lines of action gives the "center of compression". In this way, a catalogue of 
singular solutions can be built up. A number of such solutions has been given by Boussinesq 
[4] and Dougall [5], and, more recently, using other methods, Mindlin [(>] has given an extensive 
list of such "nuclei of strain", as they have been called by Love. 

By superposition of these nuclei of strain in our infinite solid, one may hope to be able to 
solve problems for domains with boundaries, and this has been done in a number of important 
cases, e. g., the half space under concentrated normal load on its boundary, or concentrated 
load on a line, or over an area such as a circle or rectangle [7]. 

The Kelvin integrals will be used here to obtain the displacements due to a concentrated 
force applied along, and normal to a circle in an infinite solid. While no mathematical features 
of particular novelty arise in the analysis, the solution does exhibit an interesting use of elliptic 
integrals. The problem has been undertaken as a first step in the application of such methods 
to problems dealing with bodies having axial symmetry; in particular, it is hoped to use the 
result obtained here as a first step in obtaining the solution due to a ring of pressure applied to 
the surface of an elastic cone. It is therefore natural to carry out the work in spherical, rather 
than the somewhat simpler cylindrical, coordinates, and this has been done. 

2. Kelvin's Integral of the Navier Equations 

Consider an elastic solid of infinite extent, and let its material points be referred to rec- 
tangular cartesian coordinates x tj i= 1,2,3. If a body force with components Fj(x ly x 2 ,x d )j 



* This research was supported by the United States Air Force, through the Office of Scientific Research of the Air Research and Develop- 
ment Command. 

365 



,7 = 1,2.3, per unit volume acts within a region V of the solid, and u j (x 1 ,x 2 ,x 3 ),j= 1,2,3 represent 
the displacement components in the solid, then the Navier equations in V are 



M V 2 ^+(\+ M )^-A: 



■F iy i=l,2,3, 



(2.1) 



where X and /x are Lame's elastic constants, and A is the sum of the extensional strains. The 
right side of (2.1) vanishes outside of V. 

Lord Kelvin [1,3] has given a solution of these equations, valid throughout the solid, which 
represents the effects of the body force, viz., 



u i{x )=AJ v [B ?f-± (1) ± (,,-^)] dV, 



i=l,2,3 



(2.2) 



where £,-, ^'=1,2,3 are variables of integration, x and £ have been used to indicate the three 
variables Xj and £ ; -, and 



3=1 



A=- r 



X+M 



B= 



X+3/x 



~8m7t(X+2ai) X+m 

3. The Kelvin Integrals in Spherical Coordinates 



(2.3) 
(2.4) 



It will be convenient for our purposes to obtain the spherical coordinate equivalents of 
(2.2), i. e., to obtain displacement components in the spherical coordinate directions expressed 
in terms of spherical coordinates. 

Let r,6, and <£ be spherical coordinates related to the x t coordinates by 

x x =r sin 6 cos </>, 

x 2 = r sin 6 sin <£, (3.1) 

x s = r cos 0, 
as indicated in figure 1. If Vi(x) denote the components of a vector in the x t coordinate system 




Figure 1. Concentrated force applied to a ring. 
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and V r (rj9j<t>), Ve(r,d,(f)), V^rfirf) its components in spherical coordinates, then it is easily 
shown either by tensor methods or by elementary geometry that 

V r =(V 1 cos <t>+V 2 sin 0) sin d+V 3 cos 0, 

Ve=(V x cos <t>+V 2 sin 0) cos 0— V z sin 0, (3.2) 

T^ = — Fi sin 0+ F 2 cos 0. 

Inversion of this system gives 

Vi = (V r sin d-j-Ve cos 0) cos <i>—V<t> sin 0, 

T7 2 =(V r r sin 0+ Ve cos (9) sin + T^ cos 0, (3.3) 

V z = V r cos 0— V* sin 0. 

If we let u r , uo, u<f, be the components of the elastic displacement in spherical coordinates, 
then they are connected with the cartesian components u t through (3.2) with V r , etc., replaced 
by u r , etc., and V t replaced by u t . If we suppose this done, and the expressions (2.2) for u t 
inserted in the result, then we obtain, after some simplification, 

v r (rfi 9 4>)=AJSB ^+± [rF r - P F r ][r-p cos co] } dV, 

Ve(r,6,<i>)=A ( S B -w+~L [rF r - P F r ] [sin cos X-cos sin X cos (0-/x)] V dV, (3.4) 

^(r,0,0)=^J /# §+|a [rF r -p^] sin X sin (0- M ) \ dV, 

where p, X, /x are the values of r, 0, 0, respectively, associated with the integration variables £*. 
It is easy to verify that 

jB 2 =r 2 + P 2 -2rp cos co, (3.5) 

cos co = cos 9 cos X+sin sin X cos (0 — //), (3.6) 

i. e., R is the distance between the fixed point (r, 9, 0) at which the displacements are to be 
calculated and the variable point of integration (p, X, p), while co is the angle between the 
radii vectors r and p. We let i^ r , etc., be the components_of the body force in spherical coor- 
dinates. The arguments of F r in (3.4) are p, X, p, but F r etc., are used to designate certain 
terms which are as yet incompletely expressed in spherical coordinates, viz, 

F T =[Fi(S) cos 4>+F 2 (£) sin 0] sin 6+F z (£) cos 0, 

^=[^(f) cos ct>+F 2 (0 sin 0] cos 0-^ 3 (f) sin 0, 

F^-F.iO sin 0+^ 2 (?) cos 0. 

But the jFi(£) are given in terms of F r (p, X, jjl), etc., by (3.3) with the obvious replacements, 
and hence we find that 

F r =F r cos co— Fett+F^ sin 9 sin (0— /z), 

Fe=-F T K + F e t E- ] rF 4> cos sin (0-/z), (3.7) 

F$= — (F r sin X4-.F0 cos X) sin (</>— p) J rF 4> cos (0— /z), 
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where the arguments of F r , etc., are p, X, /x, and 

E=sin 6 sin X+cos 6 cos X cos (/x— </>), 

A = sin 6 cos X— cos sin X cos (xx— <£), (3.8) 

S2=cos 6 sin X— sin cos X cos (/x— <£). 

With the insertion of (3.7) in (3.4), we obtain the displacements in spherical coordinates 
obtained from the Kelvin integrals: 

v r =A J ^ -^ [(F r cos w—FoQ+Ft sin (9 sin (</>— m)] 

+ r-pwa co ^ cog co _ p) _^ rl2+r ^ sin ^ gin (0 _ M )]\ ^ F? 

A Sv\ i [~AF r +SF,+^ cosjTsin (</>-,*)] +^ [F r (r cos co-p) 

--FVrQ+rjFV sin B sin («—/*)] J*dF, (3.9) 



Vd 



v*=--A j ^| -^ [—(F r sin X+^0 cos X) sin (t—ri+F* cos (<£— /x)] 



— -^3 sin X sin (0— /x) [F r (r cos co— p) —Fertt+rFj, sin sin (0— /x) r^F. 

4. The Displacement Integrals for a Ring of Concentrated Force 

The concept of a concentrated force F {c) acting at a point P of an infinite body may be 

arrived at in the following way. We first define a distributed body force which vanishes out- 

— ^ 

side some volume AV containing P in its interior, and inside AT 7 is some continuous F which 
has the property that 



? w = f 
Ja 



FdV. 

AV 



Now AV is let tend to zero, always with P in its interior, and F allowed to increase in such 
a way that the value of the integral maintains the constant value F ic) . Thus 

^ (c) (P)=lim f FdV. (4.1) 

AV-^OJAV 

We wish to consider a continuous distribution of concentrated force along a circle (C) 
whose center lies on the x s axis, and whose plane is parallel to that of the Xi and x 2 axes. (See 
fig. 1.) In spherical coordinates we take the equation of the circle to be 

r=a, d=a, 

where 0<a<oo j 0<a<7r. Suppose F{<j>) = {F T {<t>), F d (<f>),0) is a continuous distribution of 
force defined in the region a— Ar<><^+Ar, a— A6<^9<^a-r AO M~~ A$<^<O + A0 enclosing the 
point P: (a,a,/x). For simplicity, let us consider only the component Fe. In general, because 
the ^-direction varies from point to point in AV, this distribution of force will have a nonzero 
resultant force and moment with respect to the coordinate directions at (<z,a,xx). It is easy to 
show, however, that if the resultant Fe {c) in the ^-direction is held fixed as we let Ar, A0, A<j> tend 
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to zero, then we obtain in the limit only a force F e (c) acting in the ^-direction. The same remark 
holds for the component F T (<$>). In this way we obtain a concentrated force7 /(c) (jd) = (F r ic \ (p), 
F e {c) (m)i 0), acting at 4>=p on the circle r=a, 6=a. This force lies in the plane </>=m, and is 
therefore normal to the ring. 

The displacements corresponding to the distribution F=(F r ,Fe,0) in AV are seen from 
(3.9) to be 

v r =A < -7) [F r cos co— F e ti]+j^ [r—p cos co] [F r (r cos co — p) — ri}F e ] VdV, 
V6=A Lv{j{ [~ A ^ + S ^ +jr,Fr(r cos u-p)-rVF e ]\dV, (4.2) 

V4>=A J -j — j y [F r sin \+F e cos X]+"4 sin \[F r (r cos o)—p)—rUF e ] V sin (0— /x) r/F, 
where (ZT^p 2 sin X <7p r/X dp. Each of these expressions consists of a sum of terms of the form 

where F(4>) is either F r (<$>) or F e (<t>), and / represents its coefficient. It is not difficult to show 
that if (r,0,</>) =^(a,a,p), then under the assumed behavior of F in (4.1), 

lim f /(rA*;p,V)WrfH(^M;M,^ 1(c) (^ (4.3) 

af->o Jav 

Accordingly, if we now consider a distribution of force of strength F {c) per unit of arc length 
along the circle {C), we obtain for the displacements at (r } 6,<f>) the expressions 

v r =--A pTF\ B l F r ic) (t*) cos a>*-MV c >( M )] + 

-p-2 [r— a cos co*] [7V c) (m)(^ cos co*— a)~rft*7V c) 0u)] f d/*, 

'"^^Jo 2T i{ B[ ~^ (C)(M)A * + ^ (C)(;u)H * 1+ (4.4) 

^[F r «G*)(r cos co*-a)-rl2*F/ c HM)]| ^ 

t?^=-4 f 2 'J- /-7i[7 1 r (c) (/x) sin a+7V c) ( M ) cos <*] + 

5L™2f[f r c)( M )( rC08 w ,- a )-rO»F# w (M)]} sin (*-/*)<//*, 



where * indicates evaluation of the quantity to which it is attached at p=a, X=a, i. e., 

B*=r 2 -\-a 2 —2ar cos co*, 
cos co*=cos cos a + sin sin a cos (/*— </>), 

E*=sin sin a+cos cos a cos (/*— </>), (4.5) 

A*=:sin cos a— cos sin a cos (yu— </>), and 
S2*=cos sin a — sin cos a cos (/*— </>). 

The evaluation of the integrals occurring in (4.4) will be difficult, or impossible, in finite 
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terms, unless F 7 and Fq have the form of trigonometric polynomials. More general functions 
F r and Fe may be replaced by their Fourier series expansions. This leads to an infinite series 
for each displacement component, the terms of which contain integrals e valuable in terms of 
complete elliptic integrals of the first and second kinds. The process of such evaluation 
involves the use of recursion formulas which lead to extremely long and cumbersome expres- 
sions. We shall therefore restrict our consideration to the case in which the components 
F r ic) and F e {c are constant. This does not mean that the direction of the force is the same at every 
point of the ring, but that at every point it lies in the plane of the point and the x 3 axis, and has 
the same orientation in this plane. (See fig. 1.) If we denote the magnitude of this force by 
P and its angle with the r-direction by /3, where the positive sense of (3 is the same as that of 6, 
then F r =P cos /3 and_]7^=P sin 0, and (4.4) may be written as 



v r =AP 



V6 



-£-1 Bgi+-£2 (f~ a cos o*)(rgi—a cos w*) Vdn, 
=AP I -jjT 1 Bg 2 +-^ k*(rg x —a cos co*) Vdfi, 



(4.6) 



where 



<7i(0,0,a,jS,/z) = cos (a+/3) cos0+sin (a+/3) sinflcos (/*—<£), 
g 2 (d,<t>,a,(3,/u) = —cos (a+/3) sin0+sin (a+/3) cosflcos (/z -—<£). 



(4.7) 



The third integral of (4.4) vanishes, since the integrand is an odd function with respect to 
li=ir. It is now convenient to introduce nondimensional quantities 77 and h defined by 



Thus 

and (4.6) become 



r=r)d, B*=ah. (4.8) 

h 2 =l+r ] 2 -2r 1 cosu*, (4.9) 

^r=— 1 jA Bg x +^ (y— cos w*)(??#i— cos co*) >dn, 

(4.10) 
ve=— J ^< Bg 2 +j^aA*(r)g 1 — cos «*) Vd». 

5. Evaluation of the Integrals 

When the expressions (4.5) and (4.7) are inserted for A*, g x and g 2 , it will be seen that 
v r and ve consist of sums of terms of the form 



,/o [l+fl*-2l|COS«*]» 


(5.1) 


with coefficients which are functions of a, 13, 6, and »j. If we let 


( 


P 2 =^Lk=1 + >7 2 -2>7 cos (a+6), tf=hl m =l + V *-2r, COS (a-0), 


(5.2) 


then 




W=tf sin 2 M T^+2 2 cos 2 M T*, 
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and it is easy to show that 

■l™ — 2 I p - T=jj' (5.3) 

Jo p 2 sin 2 |+? 2 cos 2 | p 

These integrals are evaluable in terms of complete elliptic integrals of the first and second 
kinds [8] ; specifically, 

^[^(j*)-Mi*)} 
^K£-Xl'*)-I^)} 

where & = [(£> 2 — # 2 )/p 2 ]^ is the modulus of the functions Z? and F, and fc 2 +£' 2 =l. Using t hese 
quantities, the final form of the displacements is found to be 



JLAJP f 



B{K x F-2 sin sin (a+0)) ^ 



p a 



■- 2 Ijj-jKrf+lcos &—y\K^K 2 ^P- sin 2 6 sin a sin (a+/3) J 

— p-( s i n [ sm « cos £+?? 2 sin («+£)] — 77 sin d[K\ sin a+2JL 2 sin (a+jS)] ) \> > 
Ve=^ffJ,B(K 3 F-2 cos 6 sin («+0) p)+^[^ 2 ([('os/3-^]if 4 ~^ sin 6 cos sin (a+0)) 

+T2" (~~ cos sin a cos P+rj[Ki cos sin a-\-K A sin sin (a+/3)]) V 

where we have introduced the abbreviations 

1+fc' 2 
i^!=cos cos (a+$)-\ ~ — sin sin (a+0), 

K 2 =cos cos a-\ — jtj — sin sin a, (5.6) 

7£ 3 =— sin cos (a+/3)H — p — cos sin (a+/3), 
1+&' 2 

i£ 4 = — sin cos a+— p — cos sin a. 

In obtaining the results (5.5) the ring has been taken in the most general position in which 
the force specified acting on it can be written in terms of one coordinate. If we place the ring 
in a more special position, viz, in the x x x 2 plane, center at the origin, then (5.5) and (5.6) be- 
come somewhat simpler. Thus 
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Vr= ^{ B i K ' F ~ 2 sin 6 cos ^ I)+^[f^ iV+[C0S p-^'^- 


4&' 2 


- sin 2 cos (3 j 


2F 

— - p- (sin cos 0[1+*? 2 ]— 77 sin [If J +1^2 cos 


•"]}■ 




♦ 


%=4 ^{ B (^^- 2 cos e cos ^ i) + l [jS( [cos "-^a*; 


4fc' 2 


) sin 


cos cos j8 J 


22^ 

+-p (—cos cos jft+^iTJ cos 0+i?i sin cos 


«,]} 






where 

K[=— cos sin /H p — sin cos ($, 








i^=^sinft 






(5.8) 


1+&' 2 
i^3=sin sin /H — p — cos cos /3, 








14-&' 2 








6. A Limiting Form of the Displacements 






Let the ring be located as above, and let /3=7r/2. The fore* 
is now normal to the x x x 2 plane and has the sense of decreasing x : 


3 P at each point of the ring 
. The distribution thus has 


a resultant force 2waP=P. At distances far removed from the ring, it ] 


may 


be expected that 


the distribution on the ring will have essentially the same effec 
located at the origin. We show that this is the case. 

If we let the point (r,0,<£) recede to infinity, then lim 1/97 = 0, 


t as a concentrated force P 
and we see from (5.2) that 


for r^a 






(6.1) 


where lim 61= lim e 2 =0. Moreover, 








1A p 2 —q 2 4 sin 
V v 






(6.2) 


from which we see that lim k=0, and hence that for r^>a 






j 


<!'*)=!+-< F &)~2+«> 






(6.3) 


pk 2 =4 sin 0+e 5 






(6.4) 


where lim €* = 0, i=3, 4, 5. 

r-»co 








Using these results, it can be shown that for r^>a 








a ,t> . ^ \ 7=; cos 8 A -ryTi sin 

v r ^—A(B+l)P , vec^ABP 

v r 






(6.5) 
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If these are converted to cartesian coordinates, we find that for 0>a 

v^-AP x -f, v^-AP^ V ^-AP(B\+?). (6.6) 

The right-hand members are the displacements produced in an infinite elastic medium by a 
concentrated force P at the origin directed along the negative x z axis [3]. 
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